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It has long been thought that the sensitivity of laser interferometric gravitational-wave 
detectors is limited by the free-mass standard quantum limit, unless radical redesigns of the 
interferometers or modifications of their input/output optics are introduced. Within a fully 
quantum-mechanical approach we show that in a second-generation interferometer composed 
of arm cavities and a signal recycling cavity, e.g., the LIGO-II configuration, (i) quantum 
shot noise and quantum radiation-pressure-fluctuation noise are dynamically correlated, (ii) 
the noise curve exhibits two resonant dips, (iii) the Standard Quantum Limit can be beaten 
by a factor of 2, over a frequency range A/// ~ 1, but at the price of increasing noise at 
lower frequencies. 
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I. INTRODUCTION 



Several laser interferometric gravitational- wave (GW) detectors |jj (interferometers for short), sensitive to 
the high-frequency band 10 — 10 3 Hz, will become operative within about one year. In the first generation of 
these interferometers the Laser Interferometer Gravitational Wave Observatory (LIGO), TAMA and Virgo 
configurations [] are characterized by kilometer-scale arm cavities with four mirror-endowed test masses, 
suspended from seismic-isolation stacks. Laser interferometry is used to monitor the relative change in the 
positions of the mirrors induced by the gravitational waves. The Heisenberg uncertainty principle, applied to 
the test masses of GW interferometers states that, if the relative positions are measured with high precision, 



then the test-mass momenta will be perturbed. As time passes, the momentum perturbations will produce 
position uncertainties, which might mask the tiny displacements produced by gravitational waves. If the 
momentum perturbations and measurement errors are not correlated, a detailed analysis of the above process 
gives rise to the standard quantum limit (SQL) for interferometers: a limiting (single-sided) noise spectral 
density S^ h = 8h/(mQ 2 L 2 ) for the dimensionless gravitational-wave signal h(t) — AL/L Here m is 
the mass of each identical test mass, L is the length of the interferometer's arms, AL is the time evolving 
difference in the arm lengths, is the GW angular frequency, and h is Planck's constant. 



1 GEO's optical configuration differs from that of LIGO/TAMA/Virgo — it does not have Fabry- Perot cavities in 
its two Michelson arms, and the analysis made in this paper does not directly apply to it. However, we note that 
GEO, already in its first implementation, does use the 'signal recycling' optical configuration with which this paper 
deals. 
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The concept of SQL's for high-precision measurements was first formulated by Braginsky ||. He also 
demonstrated that it is possible to circumvent SQL's by changing the designs of the instruments, so they 
measure quantities which are not affected by the uncertainty principle by virtue of commuting with them- 
selves at different times [|| , j| - as for example in speed- meter interferometers || , which measure test- mass 
momenta instead of positions. Interferometers that circumvent the SQL are called quantum-nondemolition 
(QND) interferometers. Since the early 1970s, it has been thought that to beat the SQL for GW interfer- 
ometers the redesign must be major. Examples are (i) speed-meter designs |^| with their radically modified 
optical topology, (ii) the proposal to inject squeezed vacuum into an interferometer's dark port ||, and (iii) 
the proposal to introduce two kilometer-scale filter cavities into the interferometer's output port |?| so as to 
implement frequency-dependent homodyne detection j§] . Both (ii) and (iii) intend to take advantage of the 
non-classical correlations of the optical fields. These radical redesigns require high laser power circulating in 
the arm cavities (> 1 MW) and/or are strongly susceptible to optical losses which tend to destroy quantum 
correlations. In order to tackle these two important issues, Braginsky, Khalili and colleagues have recently 
proposed the GW "optical bar" scheme J^|, where the test mass is effectively an oscillator, whose restoring 
force is provided by in-cavity optical fields. For "optical bar" detectors the free-mass SQL is no more relevant 
and one can beat the SQL using classical techniques of position monitoring. Moreover, this scheme has two 
major advantages: It requires much lower laser power circulating in the cavities [||, and is less susceptible 
to optical losses. 

Research has also been carried out using successive independent monitors of free-mass positions. Yuen, 
Caves and Ozawa discussed and disputed about the applicability and the beating of the SQL within such 
models 10 . Specifically, Yuen and Ozawa conceived ways to beat the SQL by taking advantage of the 
so-called contractive states jH)]]. However, the class of interaction Hamiltonians given by Ozawa arc not 
likely to be applicable to G W interferometers (for further details see Ref . JTT] ) . 

Recently, we showed Jl^] that it is possible to circumvent the SQL for LIGO-II-type signal-recycling (SR) 
interferometers [[l3 14 1. With their currently planned design, LIGO-II interferometers can beat the SQL by 
modest amounts, roughly a factor two over a bandwidth A/^/. ^| It is quite interesting to notice that the 
beating of the SQL in SR interferometers has a similar origin as in "optical bar" GW detectors mentioned 
above ||. 

Braginsky and colleagues fl6|] ; building on earlier work of Braginsky and Khalili [Q, have shown that 
for LIGO-type GW interferometers, the test-mass initial quantum state only affects frequencies < 1 Hz, 
the dependence on the initial quantum state can be removed filtering the output data at low frequency. 
Therefore, the SQL in GW interferometers is enforced only by the light's quantum noise, not directly by the 
test mass. As we discussed in |l^] , and we shall explicitly show below, we can decompose the optical noise of a 
SR interferometer into shot noise and radiation-pressure noise, using the fact that they transform differently 
under rescaling of the mirror mass m and the light power I a . As long as there are no correlations between the 



2 If all sources of thermal noise can also be pushed below the SQL. The thermal noise is a tough problem and for 
current LIGO-II designs with 30 kg sapphire mirrors, estimates place its dominant, thermoelastic component slightly 
above the SQL 
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light's shot noise and its radiation-pressure-fluctuation noise, the light firmly enforces the SQL. This is the 
case for conventional interferometers, i.e. for interferometers that have no SR mirror at the output dark port 
and a simple homodyne detection is performed (the type of interferometer used in LIGO-I/TAMA/Virgo). 
However, the SR mirror Jl3|,|l4| (which is being planned for LIGOTlf] as a tool to reshape the noise curve, 
and thereby improve the sensitivity to specific GW sources (ItJ ) produces dynamical shot-noise / back- 
action-noise correlations, and these correlations break the light's ability to enforce the SQL. These dynamical 
correlations come naturally from the nontrivial coupling between the test mass and the signal-recycled optical 
fields, which makes the dynamical properties of the entire optical-mechanical system rather different from 
the naive picture of a free mass buffeted by Poissonian radiation pressure. As a result, the SQL for a free test 
mass has no relevance for a SR interferometer. Its only remaining role is as a reminder of the regime where 
back-action noise is comparable to the shot noise. The remainder of this paper is devoted to explaining these 
claims in great detail. To facilitate the reading we have put our discussion of the dynamical system formed 
by the optical fields and the mirrors into a separate, companion paper |11|. 

The outline of this paper is as follows. In Sec. [n] we derive the input-output relations for the whole optical 
system composed of arm cavities and a SR cavity, pointing out the existence of dynamical instabilities, 
and briefly commenting on the possibility and consequences of introducing a control system to suppress 
them. In Sec. Ill we evaluate the spectral density of the quantum noise. More specifically, in Sec. Ill A we 
discuss the general case, showing that LIGO-II can beat the SQL when dynamical correlations between shot 



noise and radiation- pressure noise are produced by the SR mirror. In Sec. IIIB, making links to previous 
investigations, we decompose our expression for the optical noise into shot noise and radiation- pressure noise 
and express the dynamical correlations between the two noises in terms of physical parameters characterizing 



the SR interferometer; in Sec. HI C we specialize to two cases, the extreme signal-recycling (ESR) and 
extreme resonant-sideband-extraction (ERSE) configurations, where dynamical correlations are absent and 
a semiclassical approach can be applied jl3|,[l4|. In Sec. IV we investigate the structure of resonances of 
the optical-mechanical system and discuss their link to the minima present in the noise curves. Finally, 
Sec. deals with the effects of optical losses, while Sec. ^ summarizes our main conclusions. The Appendix 
discusses the validity of the two-photon formalism in our context. 



II. SIGNAL-RECYCLING INTERFEROMETER: INPUT-OUTPUT RELATIONS 

In Fig. |l| we sketch the SR configuration of LIGO-II interferometers. The optical topology inside the 
dashed box is that of conventional interferometers such as LIGO-I/TAMA/Virgo, which are Michelson 
interferometers with Fabry-Perot (FP) arm cavities. The principal noise input and the signal and noise 
output for the conventional topology are a and di in Fig. |l|. In a recent paper, KLMTV have derived 
the input-output (ci — di) relations for a conventional interferometer at the output dark port, immediately 



3 The LIGO-II configuration will also use a power-recycling cavity to increase the light power at the beamsplitter. 
The presence of this extra cavity will not affect the quantum noise in the dark-port output. For this reason we do 
not take it into account. 



3 




signal recycling 
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FIG. 1. Schematic view of a LIGO-II signal-recycling interferometer. The interior of the dashed box refers to the 
conventional interferometer; a and di are the input and output fields at the beam splitter's dark port; a,i and bi are 
the vacuum input and signal output of the whole optical system. The laser light enters the bright port of the beam 
splitter. The arrows close to arm cavities' extremities indicate gravitational-wave-induced mirror displacements. 

after the beam splitter, within a full quantum mechanical approach. In this section we shall derive the 
input-output (dj — bi) relations for the whole optical system at the output port, i.e. immediately after the 
SR mirror, and shall evaluate the corresponding noise spectral density. 

As we shall see, a naive application of the Fourier-based formalism developed in (t]] gives ill-defined input- 
output relations, due to the presence of optical-mechanical instabilities. These instabilities have an origin 
similar to the dynamical instability of a detuned FP cavity induced by the radiation-pressure force acting 
on the mirrors, which has long been investigated in the literature ]l8|-p0|. To suppress the growing modes 
and make the KLMTV's formalism valid for SR interferometers, an appropriate control system should be 
introduced. The analysis of the resulting interferometer plus controller requires a detailed description of the 
dynamics of the whole system and for this we have found Braginsky and Khalili's theory of linear quantum 
measurement |I| very powerful and intuitive. We analyzed the details of the dynamics in an accompanying 
paper p!l| , showing in particular that the results derived in this section by Fourier techniques, notably the 
noise spectral density curves, are correct and rigorously justified. 



A. Naive extension of KLMTV's results to SR interferometers 

As in Ref. we shall describe the interferometer's light by the electric field evaluated on the optic axis 
(center of light beam) and at specific, fixed locations along the optic axis. Correspondingly, the electric fields 
that we write down will be functions of time only: all dependence on spatial position will be suppressed from 
our formulae. 

The input field at the bright port of the beam splitter, which is assumed to be infinitesimally thin, is a 
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Quantity 



Symbol Sz value for LIGO-II 



Quantity 



Symbol & value for LIGO-II 



Light power at beam splitter 
SQL for GW detection 
Laser angular frequency 
End-mirror mass 
SR cavity length 
SR mirror transmissivity 
Arm-cavity power loss 
Photodctcctor loss 



; S^ QL - 4 X 10" 48 /Hz 
= 1.8 x 10 15 sec" 1 
m — 30 kg 
I 10 m 
r (amplitude) 
e = 0.01 
A PD = 0.1 



Light power to reach SQL 
Arm-cavity half bandwidth 
GW angular frequency 
Arm-cavity length 
Internal arm-cavity mirror transmissivity 
SR cavity detuning 
SR power loss 



/sql - 1.0 x 10 4 W 
7 = Tc/4L = 2ir x 100 sec" 
Q 

h = 4 km 
T = 0.033 (power) 
<f> 

A S r = 0.02 



TABLE I. Summary of LIGO-II parameters l2B 



carrier field, described by a coherent state with power I Q and angular frequency uj . We assume that the 
arm-cavity end mirrors oscillate around an equilibrium position that is on resonance with the carrier light. 
This means that there is no zeroth-order arm-cavity detuning (see the paper of Pai et al. [ p0| for a critical 
discussion of this assumption). Our most used interferometer parameters are given in Table |l] together with 
the values anticipated for LIGO-II. 

We denote by /gw = ^/27r the GW frequency, which lies in the range 10— 1000 Hz. Then the interaction 
of a gravitational wave with the optical system produces side-band frequencies w ± J7 in the electromagnetic 
field at the output dark port. For this reason, similarly to KLMTV 0, we find it convenient to describe the 
quantum optics inside the interferometer using the two-photon formalism developed by Caves and Schumaker 
p2| , p3| . In this formalism, instead of using the usual annihilation and creation operators for photons at 
frequency to, we expand the field operators in terms of quadrature operators which can simultaneously 
annihilate a photon at frequency u> = lo + CI while creating a photon at frequency lo = lo — Q (or vice 
versa) . 

More specifically, the quantized electromagnetic field in the Heisenberg picture evaluated at some fixed 
point on the optic axis, and restricted to the component propagating in one of the two directions along the 
axis is: 

E(t) = jf" ^ e+l "] ^ ■ 

Here A is the effective cross sectional area of the laser beam and c is the speed of light. The annihilation 
and creation operators a u , in Eq. ( |2.l| ), which in the Heisenberg picture are fixed in time, satisfy the 
usual commutation relations 

\a u ,a u .]=0, [o£,o£,]=0, [a u , a*,] = 2n5(u - u') . (2.2) 

Henceforth, to ease the notation we shall omit the hats on quantum operators. Defining the new operators 
(see Sec. IV of Ref. §§ g) 



Our notations are not exactly the same as those of Caves and Schumaker |22|,|23|, the correspondence is the 
following (ours — > Caves- Schumaker): u>o — » fi, fi — » e, eiu; ±n — > 0±, a± — > X±a±, a±.2 — » Qfi,2- We refer to Sec. IV 
B of Q for further details. 
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and using the commutation relations (2.2), we find 

[a + ,a\_,] = 2n5{n-n') + [o_, 4,] = 2-kS(Q - Si') (l - ^) , (2.4) 

[a + , a + i] = = [a_, a_/] , [a\_, a',,] = = [at, 4/] , [a + , a_/] = = [a + , 4/] , (2-5) 

where a±< stands for a±(fl'). Because the carrier frequency is lu ~ f0 15 s _1 and we are interested in 

frequencies H/2tt in the range 10 10 3 Hz, we shall disregard in Eq. (2J3) the term proportional to 

Q/u . (In the Appendix we shall give a more complete justification of this by evaluating the effect the term 
proportional to fl/to would have on the final noise spectral density.) We can then rewrite the electric field, 
Eq. (O), as 



£'(/) : : J^±e~^ 4 [ +00 (a + (ri)e-^ + a_(n)e^)^ + h.c, (2.6) 



Ac Jo 2vr 

where "h.c." means Hermitian conjugate. Following the Caves-Schumaker two-photon formalism J2^,^3|, we 
introduce the amplitudes of the two-photon modes as 

«+ + 4 a+~ a - , n 

ai and a 2 are called quadrature fields and they satisfy the commutation relations 

[ai, 4'] = -[02, 4'] = 27ri(5(n - Q!) , 

[ai,a\,] = = [ai,a v ] , [a 2 , 4'] = = [«2, ay] ■ (2.8) 



Expressing the electric field ( |2.6| ) in terms of the quadratures we finally get 

E(a,i;i) = cos(uj t) E^a^t) + sin(uj t) E 2 (a 2 ;t) , (2.9) 

with 

Ej (a r ,t) = ^ + °°(a, e"™ + a] e^) f£ , = 1,2. (2.10) 

Note (as is discussed at length by BGKMTV (l6) and was previewed by KLMTV (footnote 1 of Ref. 0)), 
that, Ei(t) and E2(t) commute with themselves at any two times t and t' , i.e. [Ej(t) t Ej(t')] — 0, while 
[Ei(t) t Ezit')] ~ i<5(£ — £'). Hence, the quadrature fields Ej(t) with j = 1,2 are quantum-nondemolition 
quantities which can be measured with indefinite accuracy over time, i.e. measurements made at different 
times can be stored as independent bits of data in a classical storage medium without being affected by 
mutually- induced noise, while it is not possible to do this for Ei(t) and E%(t) simultaneously. As BGKMTV 
fl6f emphasized (following earlier work by Braginsky and Khalili Q), this means that we can regard Ei(t) 
and E 2 (t) separately as classical variables - though in each other's presence they behave nonclassically. 

For GW interferometers the full input electric field at the dark port is E(ci\t) where c\ and c 2 are the 
two input quadratures, while the output field at the dark port is E{di]t) 1 with d\ and d 2 the two output 
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quadratures (see Fig. [j]). Assuming that the classical laser- light input field at the beam splitter's bright 
port is contained only in the first quadrature]^ and evaluating the back-action force acting on the arm-cavity 
mirrors disregarding the motion of the mirrors during the light round-trip time (quasi-static approximation) , 
KLMTV derived the following input-output relations at side-band (GW) angular frequency Q 

d 1 = c 1 e 2i l 3 , d 2 = (c 2 - ICa) e 2tfl + V2IC e lf3 , (2.11) 

"SQL 

where 2/3 = 2 arctan n/7 is the net phase gained by the sideband frequency f2 while in the arm-cavity, 
7 = Tc/AL is the half bandwidth of the arm-cavity (T is the power transmissivity of the arm-cavity input 
mirrors and L is the length of the arm cavity), h is the Fourier transform of the gravitational- wave field, 
and hsQL is the SQL for GW detection, explicitly given by 



where m is the mass of each arm-cavity mirror. The quantity JC in Eq. ( 2.11 ) is the effective coupling 
constant, which relates the motion of the test mass to the output signal, 



Finally, I a is the input light power, and Isql is the light power needed by a conventional interferometer to 
reach the SQL at a side band frequency CI = 7, that is 

(See in Table | the values of the interferometer parameters tentatively planned for LIGO-II |2^].) We shall 
now derive the new input-output (a, — 6j) relations including the SR cavity. We indicate by I the length of 
the SR cavity and we introduce two dimensionless variables: 4> = [w Z/c] mo d27r, f\ the phase gained by the 
carrier frequency oj while traveling one way in the SR cavity, and <3? = [Cll j '(^xaoA2n the additional phase 
gained by the sideband with GW frequency ft (see Fig. [j]). Note that we are assuming that the distances 
from the beam splitter to the two arm-cavity input mirrors are identical, equal to an integer multiple of the 
carrier light's wavelength, and are negligible compared to I. 

Propagating the output electric field E(di\t) up to the SR mirror, and introducing the operators e» and 
fi which describe the fields that are immediately inside the SR mirror (see Fig. 1), we get the condition 

E{f i ;t)=E(d i ;t-~\ , (2.15) 



5 For the KLMTV optical configuration and for ours, only a negligible fraction of the quantum noise entering the 

bright port emerges from the dark port. 

6 The description of a SR interferometer beyond the quasi-static approximation po|Jl9[ ] introduces nontrivial cor- 
rections to the back-action force, proportional to the power transmissivity T of the input arm-cavity mirrors. Since 
T ~ 0.033 (see Table [j) we expect a small modification of our results, but an explicit calculation is strongly required 
to quantify these effects. 

7 Note that u> l/c = 2-Krn + (f>, with m a large integer. Indeed, typically u> — 10 15 s _1 , I ~ 10m, hence lu 1/c 3> 1. 
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which, together with Eq. ( |2.9| ), provides the following equations 

fi = (di cos0 — d 2 sin</>) e l * , /a = {d\ smcj) + d 2 cos</!>) e** . (2-16) 
Proceeding in an analogous way for the input electric field E(cf,t), we derive 

e\ = (ci cos^ + c 2 sin</>) e _I * , e2 = (— Ci sin</> + c 2 cos0)e~ l * . (2-17) 



Note that each of Eqs. ( [2.16 ), ( [2.17 ) correspond to a rotation of the quadratures di, d 2 (or ci, C2) plus the 



addition of an overall phase. Finally, denoting by a, and 6^ the input and output fields of the whole system 
at the output port (see Fig. |]) we conclude that the following relations should be satisfied at the SR mirror: 

e 1 =rai+pf 1 , e 2 = ra 2 +pf 2 , (2-18) 
h = rfi-pa 1) b 2 = rf 2 -pa 2 , (2.19) 

where ±p and r are the amplitude reflectivity and transmissivity of the SR mirror, respectively. We use the 
convention that p and r are real and positive, with the reflection coefficient being +p for light coming from 
inside the cavity and — p for light coming from outside. In this section we limit ourselves to a lossless SR 
mirror; therefore the following relation holds: r 2 + p 2 = 1. 

Before giving the solution of the above equations, let us notice that the equations we derived so far for 
the quantum EM fields in the Heisenberg picture are exactly the same as those of classical EM fields. To 
deduced them it is sufficient to replace the quadrature operators by the Fourier components of the classical 
EM fields. The input-output relation we shall give below is also the same as in the classical case. In the 
latter we should assume that a fluctuating field enters the input port of the entire interferometer. More 
specifically, assuming a vacuum state in the input port, we can model the two input quadrature fields as two 
independent white noises. Then using the classical equations, we can derive the output fields which have 
the correct noise spectral densities. 

Solving the system of Eqs. ( 2.11 ), ( 2.16 )-( ^.19 ) gives the final input-output relation: 



bi \ _l_ 
bo) ~ M 



. C 2 i C22 / V a 2 J \D 2 J «sql 

where, to ease the notation, we have defined: 



(2.20) 



M = 1 + p 2 e 4 ^+*> - 2p e 2 ^+*) (^cos 20 + - sin 2^ , (2.21) 
C n = C22 = (1 + P 2 ) (cos 2<P + ^ sin 2(f> ) - 2p cos (2(f3 + $)) , (2.22) 



2 

Cxi = -t 2 (sin 2(j) + K, sin 2 0) , C 2 i = r 2 (sin 2$ - K. cos 2 <p) , (2.23) 

£>i = -(l + pe 2lil3+ ^) sin0, D 2 = -{-1 + p e 2l ^+^) cos0. (2.24) 

A straightforward calculation using € 3? and Cn C22 — C± 2 C 2 \ = |M| 2 , confirms that the quadratures 
bi,b\ satisfy the commutation relations (|1]), as they should since like and a\ they represent free fields. 
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Let us also observe that both the quadratures b\ and b 2 in Eq. ( 2.20 ) contain the gravitational- wave signal 
h and that it is not possible to put the signal into just one of the quadratures through a transformation that 
preserves the commutation relations of 61 and b 2 . Indeed, the most general transformation that preserves 
the commutation relations is of the form 

b 2 J \L 2 i L 22 J \b 2 . 



Lij e Sft , det Li 



(2.25) 



where a is an arbitrary phase. Because the Di are complex [see Eq. ( 2.24 )], it is impossible to null the h 
contribution either in b\ or b 2 . 

Henceforth, we limit our analysis to $ = 0, which corresponds to a SR cavity much shorter than the 
arm-cavities, e.g., I ~ 10m. We assume for simplicity that there is no radio-frequency (MHz) modula- 
tion/demodulation of the carrier and the signal [^lf ; instead, some frequency-independent quadrature 



b^= b\ sin Q + b 2 cos ( 
= J} [ e2l/3 ( Cl1 sin C + ^2i cos C) 01 +e 2l/3 (Ci 2 sinC + C 2 2 cosC)a 2 



2/Cre l/3 (Dx sinC + A; cosC) 



ISQL 



(2.26) 



is measured via homodyne detection ^] Before going on to evaluate the noise spectral density in the 
measured quadrature b^ , let us first comment on the results obtained in this section. 



B. Discussion of the naive result 



There is a major delicacy in the input-output relation given by Eq. ( |2.18 ). By naively transforming it 



from the frequency domain back into the time domain, we deduce that the output quadratures depend on 
the gravitational-wave held and the input optical holds both in the past and in the future. Mathematically 



this is due to the fact that the coefficient 1/M, in front of h and di (i = 1,2) in Eq. (2.2C), contains poles 
both in the lower and in the upper complex plane. This situation is a very common one in physics and 
engineering (it occurs for example in the theory of linear electronic networks [25|| and the theory of plasma 
waves [p6| ), and the cure for it is well known: in order to construct an output held that only depends on 
the past, we have to alter the integration contour in the inverse-Fourier transform, going above (with our 
convention of Fourier transform) all the poles in the complex plane. This procedure, which can be justified 
rigorously using Laplace transforms [p7f , makes the output signal inhnitely sensitive to driving forces in 
the inhnitely distant past. The reason is simple and well known in other contexts: our optical mechanical 



system possesses instabilities, which can be deduced from the homogeneous solution b^° m of Eqs. (2.11) 



(2.15) and (2.19), which has eigenfrequencies given by M = 0. Because the zeros of the equation M = 



8 It is still unclear what detection scheme (direct homodyne detection or RF modulation/demodulation ) will be 
used in LIGO-II. The decision will require a quantum-mechanical analysis of the additional noise introduced by the 
modulation/demodulation process, which will be given in a future paper Q. 
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are generically complex and may have positive imaginary parts [fTT| , we end up with homogeneous solutions 
that grow exponentially. f\ 

To quench the instabilities of a SR interferometer we have to introduce a proper control system. In 
we have given an example of such a control system, which we briefly illustrate here. Let us suppose that 
the observed output is b^ and we feed back a linear transformation of it to control the dynamics of the end 
mirrors. This operation corresponds to making the following substitution in Eq. ( |2.26| ): 

h^h + Cb (2.27) 

where C is some retarded kernel. Solving again for b^, we get: 

b <= TT [e 2l/3 (di sinC + C 21 cosC) 01 + e 2l0 {C 12 sinC + C 22 cosC) a 2 
Mc 



2Kre lP (D 1 sm£ + D 2 cosC) 



*SQL 



(2.28) 



simply replacing the M in (2.26) by Mc, which depends on C. Note that, by contrast with the uncontrolled 
output Eq. (2.20), the output field is no longer a free electric field, i.e. a quadrature field defined in half 
open space, satisfying the radiative boundary condition. This is due to the fact that part of it has been fed 
back into the arm cavities. Nevertheless, in the time domain, b^ commutes with itself at different times. In 



|L1| we have shown that there exists a well-defined C that makes Eq. ( 2.28] ) well defined in the time domain, 
getting rid of the instabilities. As a consequence, Mc has zeros only in the lower-half complex plane and we 
can neglect the homogeneous solution Mc b^ hom — because it decays exponentially in time. 

Finally, let us remember the important fact that the introduction of this kind of control system only 
changes the normalization of the output field. As a consequence, the noise spectral density is not affected. 
However, an extra noise will be present due to the electronic device that provides the control force on the 



end mirrors. K. Strain estimated that it can be kept smaller than about 10% of the quantum noise 28 



III. FEATURES OF NOISE SPECTRAL DENSITY IN SR INTERFEROMETERS 



In light of the discussion at the end of the last section, we shall use Eq. (2.28) as the starting point of our 
derivation of the noise spectral density of a (stabilized) SR interferometer. 

A. Evaluation of the noise spectral density: going below the standard quantum limit 



The noise spectral density is calculated as follows JjJ. Equation ( 2.28| ) tells us that the interferometer 



noise, expressed as an equivalent gravitational-wave Fourier component, is 

, _ h S QL 



A&c, (3.1) 



Quadrature operators with complex frequency can be defined by analytical continuations of quadrature operators 
with real frequency considered as analytical functions of Q. 
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where 

_ (Cn sin C + On cos Q Qi + (Ci 2 sin ( + C 22 cos Q a 2 ,„ 
C ~ r(£> lS inC + i3 2 cosC) ' 1 ' 

Then the (single-sided) spectral density S^f), with / = S1/27T, associated with the noise h n can be computed 
by the formula (Eq. (22) of Ref. @) 

i 27r<5(r! - n') $£(/) = (in|/i„(n) ^(n')|in) sym = ^<m|ft„(Q) ftt (fi') + fct (fi') /i„(n)|in) . (3.3) 

Here we put the superscript ( on to remind ourselves that this is the noise when the output is monitored 
at carrier phase C by homodyne detection. Assuming that the input of the whole SR interferometer is in its 
vacuum state, as is planned for LIGO-II, i.e. |in) = |0 Q ), and using 

<0ak4|0 a ) sym = X - 27T5(Q - fi')% , (3.4) 



(Eq. (25) of Ref. 0) we find that Eq. (3.3) can be recast in the simple form (note that Cy £ 3?): 

Q C _ h SQL (di sinC + C 2 i cosC) 2 + (C i2 sinC + C 22 cosC) 2 , , 



For comparison, let us recall some properties of the noise spectral density for conventional interferometers 
(for a complete discussion see Ref. M). To recover this case we have to take the limit (f> — *■ and p — > in 



the above equations or simply use Eq. (2.11) (in a conventional interferometer there are no instabilities). In 



particular, for a conventional interferometer, Eqs. ( 2.26| ) and ([O]) take the much simpler form 



65 onv = cos ({[a 2 + (tan ( - K.) a x ] e 2i0 } , h c ° av = e tfs [a 2 + (tan £ - K) a x ] , (3.6) 



and the noise spectral density reads 



hi 



S C r n V= ^L [l + (t&aC _ ]C) 2 ] , (3J) 

As has been much discussed by Matsko, Vyatchanin and Zubova [B| and by KLMTV M, and as we shall see 



in more detail in Sec. IIIB, taking as the output instead of the quadrature 6 2 in which all the signal h is 
encoded, builds up correlations between shot noise and radiation-pressure noise. We refer to correlations of 
this kind, which are introduced by the special read-out scheme, as static correlations by contrast with those 
produced by the SR mirror, which we call dynamical since they are built up dynamically, as we shall discuss 



in Sec. IV. The static correlations allow the noise curves for a conventional interferometer to go below the 
SQL when I Q = Isql, as was originally observed by Matsko, Vyatchanin and Zubova M. However, if £ 
is frequency independent as it must be when one uses conventional homodyne detection, then the SQL is 
beaten, «S^' conv < ^-sql, only over a rather narrow frequency band and only by a very modest amount. On the 
other hand, as Matsko, Vyatchanin and Zubova || showed, and one can see from Eq. ( |3.7| ), if we could make 



the homodyne detection angle £ frequency dependent, then choosing |?j C(fi) = arctan/C(f2), would remove 



10 Note that our definition of £ differs from the one used in 
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FIG. 2. Log-log plot of y/S h {Q)/y/sl QL (-y) versus Q/j for (i) the quadratures 6? (C = tt/2) and b% (£ = 0) with 
p = 0.9, = 7r/2 — 0.47 and I — /sql, (h) the SQL, (hi) a conventional interferometer with I = /sql, and (iv) the 
noise curve of LIGO-II @ one would obtain if shot-noise / radiation-pressure correlations were (naively) neglected. 
For LIGO-II, 7 = 2tt x 100 Hz (top axis) and y/S^^y] = 2 x 10~ 24 Hz _1/2 . These curves do not include seismic 
and thermal noises; for LIGO-II the latter is expected to be slightly above the SQL M . 



completely (in the absence of optical losses) the second term in the square parenthesis of Eq. (3.7), which is 
the radiation-pressure noise, leaving only the shot noise in the interferometer output, i.e. g£> conv = fogQ L /2/C. 
In order to implement frequency dependent homodyne detection, KLMTV have recently propose to place 
two 4km-long filter cavities at the interferometer dark port and follow them by conventional homodyne 
detection. This experimentally challenging proposal would allow the interferometer to beat the SQL at 
frequency / = 100 Hz by a factor y/S^/\Jsl QL ~ 0.24, over a band of A/-/, at light power I a = J SQL , 

and by \J S™ m / \J S^ QL ~ 0.18 if I Q ~ 3.2/sql- In conclusion, already in conventional interferometers it is 
possible to beat the SQL provided that we measure and build up proper static correlations between shot 
noise and radiation-pressure noise. 

Let us now go back to SR interferometers. They have the interesting property of building up dynamically 
the correlations between shot noise and radiation-pressure noise, thanks to the SR mirror. Indeed, even if we 
restrict ourselves to the noise curves associated with the two quadratures 6^ and b%, i.e. we do not measure 
the SR interferometer can still go below the SQL. Moreover, if the SR interferometer works at the SQL 
power, i.e. I D = Isql, as is tentatively planned for LIGO-II, then the noise curves [Eq. ( |3.5| )] can exhibit 
one or two resonant dips whose depths increase and widths decrease as the SR-mirror's reflectivity is raised. 
(We postpone the discussion of this interesting feature to Sec. 0.) These resonances allow us to reshape the 
noise curves and beat the SQL by much larger amounts than in a conventional interferometer with static 
correlations introduced by frequency-independent homodyne detection. 

More specifically, the noise spectral density, Eq. ([Te|), depends on the physical parameters which char- 
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FIG. 3. Log-log plot of ^/S h (n)/5^ L ( 7 ) versus Q/7 for the following choices of the frequency independent 
homodyne phase: £ = 0, C = C = "Y3 and £ = tt/2, with p = 0.9, = 7r/2 — 0.47 and I = /sql- The 

plot also shows the noise curve for a conventional interferometer and the SQL line. For LIGO-II, 7 = 2tv x 100 Hz 
(top axis) and y/sf^y) = 2 x 10" 24 Hz" 1 / 2 . 



acterize the SR interferometer (see Table |): the light power I a , the SR detuning 0, the reflectivity of the 
SR mirror p and the homodyne phase £. To give an example of LIGO-II noise curves, in Fig. |^ we plot the 
yj Shift) for the two quadratures (( = vr/2) and b% (( = 0), for: p = 0.9, <j) = k/2 - 0.47 and I Q = Isql- 
Also shown for comparison are the SQL line, the noise curve one would obtain if one ignored the correla- 
tions between the shot noise and radiation-pressure noise pl| , P] and for a conventional interferometer with 
lo = Isql and C = 0, explicitly given by 0: 



c*C— 0,conv 



-fSQL 



K 



(3.8) 



The sensitivity curves for the two quadratures go substantially below the SQL and show two interesting 
resonant valleys. In Fig. |^ we plot the noise curves \J Shift) f° r different values of the frequency independent 
homodyne angle (, choosing the same parameters used in Fig. g, i.e. p — 0.9, (f> — ir/2 — 0.47 and I — Isql- 
Note that the location of the resonant dips does not depend much on the angle £. This property is confirmed 



analytically in Sec. IV in the case of a highly-reflecting SR mirror, by an analysis that elucidate the underlying 



11 Before the research reported in this paper, the LIGO community computed the noise curves for SR interfer- 
ometers by (i) evaluating the shot noise S^ ot , (ii) then (naively assuming no correlations between shot noise and 
radiation-pressure noise) using the uncertainty principle S^ hot S^ p > (S^ 1 ") /4, with the equality sign to evaluate 
the radiation-pressure noise Sj^ F , (iii) then adding the two. This procedure gave the noise curve labeled "correlations 
neglected" in Fig. |; see Fig. 2 of Ref. 
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physics. 

Before ending this section, let us give an idea of the performances achievable in a SR interferometer if its 
thermal noise can be made negligible p5|| . We have estimated the signal-to-noise ratio for inspiraling binaries, 
which are among the most promising sources for the detection of GW with earth-based interferometers. The 
square of the signal-to- noise ratio for a binary system made of black holes and/or neutron stars is given by: 

Using the Newtonian quadrupole approximation, for which the waveform's Fourier transform is |/i(/)| 2 oc 
J~ 7 / 3 , and introducing in the above integral a lower cutoff due to seismic noise at fl s = O.I7 (f s ~ 10 Hz), 
we get for the parameters used in Fig. ||: 

(S/Nh ~1.83, ,j S M> *1.98, (3.10) 



(S/N) conv ■ ' (S/N) c 

where (S/N)i, (S/N)2 and (S/N) conv use for the noise spectral density either that of the first quadrature b^ 
or the second quadrature b% or the conventional interferometer, respectively. A more thorough analysis of 
signal-to-noise ratio for inspiraling binaries inevitably requires the specification of the readout scheme and 
we plan to publish it elsewhere p3| . 



B. Effective shot noise and radiation- pressure noise 

In this section we shall discuss the crucial role played by shot-noise / radiation-pressure correlations that 
are present in LIGO-II's quadrature outputs ( 2.2CQ and noise spectral densities (3.5), in beating the SQL. Our 



analysis is based on the general formulation of linear quantum measurement theory developed by Braginsky 
and Khalili Q and assumes also the results obtained in |l6| , Q . 

To identify the radiation pressure and the shot noise contributions in the total optical noise, we use the 
fact that they transform differently under rescaling of the mirror mass. Indeed, it is straightforward to show 
that in the total optical noise there exist only two kinds of terms. There are terms that are invariant under 
rescaling of the mass and terms that are proportional to 1/m. Hence, quite generally we can rewrite the 
output O of the whole optical system as [Q, |ll[] 

e>(fi) = z{Q) + n xx {Q) + l h(n) , (3.11) 

where by output we mean one of the two quadratures 6j, b^ or a combination of them, e.g., b^ (modulo a 
normalization factor) and where TZ X x is the susceptibility of the antisymmetric mode of motion of the four 
mirrors Q , given by 

fta*(fi) = ts- (3-12) 



The observables Z and T in Eq. (3.11) do not depend on the mirror masses m, and satisfy the commutation 
relations (Eq. (2.19) in Rcf. 0) 

= 0= [z(Q),z\o!)], [z(Q),^(n')\ = -2nih5(n - n') . (3.13) 
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We shall refer to Z and T as the effective shot noise and effective radiation-pressure force, respectively, 
because we have shown jTT) that for a SR interferometer the real back-action force acting on the test masses 
is not proportional to the effective radiation-pressure noise, but instead is a combination of the two effective 
observables Z and T . When the shot noise and radiation-pressure noise are correlated, the real back-action 
force does not commute with itself at different times, which makes the analysis in terms of real quantities 
more complicated than in terms of the effective ones. We prefer to discuss our results in terms of real 
quantities separately |Tl[ ] , in a more formal context which uses the description of a GW interferometer as a 
linear quantum-measurement device Q]. 

The noise spectral density, written in terms of the effective operators Z and reads 

S h = ^ { s zz + 2n ™ ® l s Fz] + n l* Srr} , (3.14) 
where the (one-sided) cross spectral density of two operators is expressible, by analogy with Eq. ( |3.3| ), as 
- S AB (n)2ir6 (ft - Q') = ^ (A(Sl)B*(tf) + B\Sl')A(Q)) . (3.15) 



In Eq. (3.14) the terms containing Szz, Sff and 3? [Sj?z] should be identified as effective shot noise, back- 



action noise and a term proportional to the effective correlation between the two noises, respectively 



Relying on the commutators (3.13) between the effective field operators one can derive H, |11| the following 



uncertainty relation for the (one-sided) spectral densities and cross correlations of Z and T\ 

Szz Sjrjr - Szr Stz > h 2 ■ (3.16) 



Equation (3.16) does not, in general, impose a lower bound on the noise spectral density Eq. ( [3.14 ). However, 



in a very important type of measurement it does, namely for interferometers with uncorrelated shot noise and 
back-action noise, e.g., LIGO-I/TAMA/Virgo. In this case Szr = = S^ z |Q and inserting the vanishing 
correlations into Eqs. ( 3.15Q , ( 3.16| ), one easily finds that the noise spectral density has a lower bound which 



is given by the standard quantum limit, i.e. 

From this it follows that to beat the SQL one must create correlations between shot noise and back-action 
noise. 

Before investigating those correlations in a SR inteferometer, we shall first show how such correlations can 
be built up statically in a conventional (LIGO-I/TAMA/Virgo) interferometer by implementing frequency- 
independent homodyne detection at some angle C @, @- By identifying in the interferometer output 



(3.6) the terms independent of m as effective shot noise and those inversely proportional to m as effective 

:™ nv and 



back-action noise, we get the effective field operators Z^ onv and ^"™ nv - 



12 We have shown jllj that as a consequence of this identification the antisymmetric mode of motion of the four 
mirrors acquires an optical-mechanical rigidity and a SR interferometer responds to GW signal like an optical spring. 
This phenomenon was already observed in optical bar detectors by Braginsky's group ^). 
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T"Conv 



(fi) 



(a 2 +01 tan 4) , y 7 ^ (Si) = — — a\ 



(3.18) 



>2K ' LhsQL 

[We remind the readers that hsQL oc \ j\fm and that JC cx 1/m.] Evaluating the spectral densities of those 



operators using Eqs. (3.15) and (3.4), we obtain the following expressions for the spectral densities and their 
static correlations: 



(Si) 



(SI) 



2JC 7T 

L2 h SQL 



S%? ( (O) = tan C = (SI) . (3.19) 



By inserting these in Eq. (3.14) and optimizing the coupling constant }C, we see that the SQL can be beaten 
for any < £ < 7r/2, i.e. whenever there are nonvanishing correlations. See Refs. (Q and ||] for further 
details. 

Let us now derive the correlations between shot noise and back action noise in SR interferometers. Because 
in this case the correlations are built up dynamically by the SR mirror and are present in all quadratures, as 



an example, we limit ourselves to the two quadratures b\ and 6 2 . Identifying in Eqs. (3.1), (3.2) the effective 
shot and back-action noise terms due to their m dependences, we obtain the effective field operators Z\, Z 2 , 
T\ and T 2 

e tP Lh SQL [at (-2p cos 2/3+ (1 + p 2 ) cos 20) + a 2 (-1 + p 2 ) sin 20] csc0 



Z 1 (Sl) 
Z 2 (Sl) 



and 



V2X t (1 + e 2i P p) 

e i/3 L/i SQL [ai(l-p 2 ) sin20 + a 2 (-2p cos2/3 + (1 + p 2 ) cos 20)] sec< 
V21C 

^i(fi) = 
r 2 (si) -- 



r(-l + e 2i P p) 

h e if} V2IC [01 (1 + p 2 ) cos + a 2 (-1 + p 2 ) sin 0] 

Lh SQL t (1 + e 2 ^ p) 

he i(3 V2lC [01 (-1 + p 2 ) cos0 + a 2 (1 + p 2 ) sin0] 



(3.20) 



SQL 



.(_ 1 + e 2 l /3 p ) 



(3.21) 



Evaluating the spectral densities of the above operators through Eqs. ( 3.15 ) and (3^4) we obtain the following 
expressions: 

1 + p A + 2p 2 cos 



h 2 2K 



L 2 h 2 SQL (l-p 2 )(l + p 2 + 2pcos2/3) ' 

h 2 2)C l + p 4 -2p 2 cos0 

IJh 2 ^ (l-p 2 )(\ + P 2 -2p cos2/3) ' 



(3.22) 



and 



<SZ 2 22 (^) = 



^Mol [4(-l + P 2 ) 2 cos 2 



-2/9 cos 2/3 + (1 + p 2 ) cos 2(j>) 2 esc 2 



2/C 

1/2 ^SQL 



(1 - p 2 ) (1 + p 2 + 2p cos 2/3) 
[4(-l + p 2 ) 2 sin 2 (f>+(-2p cos 2/3+ (1 + p 2 ) cos20) 2 sec 2 . 



2/C (1 - p 2 ) (1 + p 2 - 2p cos 2/3) 

Finally, for the correlations between the shot noise and back-action noise we get: 

h [(-1 + p 2 ) 2 - 2p(l + p 2 ) cos2/3 + 4p 2 cos 20] cot. 



(3.23) 



S^ lZl (Sl) = S Zx r x (Sl) = - 
Sf 2 z 2 (^) = Sz 2 t 2 (SI) 



(l-p 2 )(l + p 2 + 2p cos 2/3) 
h [(-l + p 2 ) 2 + 2p(l + p 2 ) cos2/3-4p 2 cos 20] tan< 
(l-p 2 )(l+p 2 -2p cos 2/3) 



(3.24) 
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These correlations depend on the sideband angular frequency fl and are generically different from zero. How- 
ever, when <j) — and <fi = n/2 the correlations are zero. We shall analyze these two extreme configurations 
in the following section. 



C. Two special cases: extreme signal-recycling and resonant-sideband-extraction configurations 

In this section we discuss two extreme cases that are well known and have been much investigated in 
the literature using a semiclassical analysis 13,14], In these two cases the dynamical correlations between 



shot noise and radiation-pressure noise are zero. This has two implications: (i) the semiclassical analysis 
and predictions Jl3|,|l4| are correct (when straightforwardly complemented by radiation pressure noise) , and 
(ii) the noise curves are always above the SQL. Of course, static correlations can always be introduced by 
measuring the quadrature b^. In these two extreme cases there are no instabilities and the input-output 
relation of the SR interferometer can be obtained from the conventional noise by just rescaling the parameter 
K [Eq. ( p3j )]. 



1. Extreme signal-recycling (ESR) configuration: (f> = 

For = 0, the gravitational- wave signal appears only in the second quadrature 62 but not in the first 
quadrature b\ (see Eq. ( 2.26| ) with £ = and n/2, respectively). Defining 

K.T 2 

K = l + pi-2p cos2/3 ' (3 - 25) 
it is straightforward to deduce that the spectral density of the noise takes the simple form 

Sf R =\- (i+JCj . (3.26) 



In the left panel of Fig. | we plot J S^ SR (n)/Sl QL (-f ) versus Q/7 for different choices of the reflectivity p. 
As we vary the reflectivity of the SR mirror the minimum of the various curves is shifted along the SQL 
line, and the shape of the noise curve change a bit because both JC and (3 in Eqs. ( 3.25| ), ( 3.26| ) depend on 



frequency. Moreover, for fi/7 3> 1 and O/7 <C 1 the curves are well above the conventional interferometer 
noise. This effect becomes worse and worse as p — > 1 and is described by the formulas 

sr K m t in 2 j SQ L n >>v gg sR (») 7 4 (i±a io ^ <<;1 _ 



5 sql (7) 47 2 \i- p J i o ' 7 ' 5^ QL ( 7 ) n 4 \1- 9) ^sql' 7 

(3.27) 

The signal-to-noise ratio for inspiraling binaries is given in this case (for p = 0.9, I a = /sql) by: 

(3.28) 

W/-'* Jconv 

Hence, this LIGO-II configuration (4> = 0) is not appealing. The noise curves could be better than the ones 
for a conventional interferometer in the range ^20 — 60 Hz, depending on the value of p, but they get worse 
everywhere else, and overall, for any p the signal-to-noise ratio for inspiraling binaries is lower than in the 
case of a conventional interferometer. 
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£2/7 Clly 
FIG. 4. Log-log plot of ^/S^ SR (Q.)/S s h QL (j) versus fi/7 for the extreme signal-recycling configuration (left panel) 
<f> = with p — 0.7, p = 0.8, p = 0.9 and I ~ Isql and for the extreme resonant-sideband-extraction configuration 
(right panel) (f> — n/2 with p — 0.7, p = 0.8 and p = 0.9, with I = Isql- Also plotted for comparison are the noise 
curve for a conventional interferometer and the SQL line. For further detail on these well known configurations, see 
Refs. 



2. Extreme resonant- sideband- extraction (ERSE) configuration: (f> — n/2 



For <j) = tt/2, using Eq. ( 2.26 ) with £ = n/2, we find that only the first quadrature b\ contains the 



gravitational-wave signal. Introducing 

— K.T 2 

K = l + p2 + 2 pcos2/3' (3 - 29) 
(which depends on frequency through both K, and (3) , we easily deduce that the noise spectral density reads 

S er SE = ^_ n_ +1 A . (3.30) 



The right panel of Fig. | shows yj S% RSE (ft) / Sl QL (7) as a function of ft/7 for different values of the reflec- 
tivity p. As for the ESR configuration discussed above, when we vary the reflectivity of the SR mirror the 
minimum of the various curves moves along the SQL line. But by contrast with the ESR configuration, for 
Q/7> 1 and ft/7 1 the curves are significantly below the conventional-interferometer noise. This effect 
becomes better and better as p — > 1 and is described by the asymptotic limits 

Sf^) 1ft 2 /1-pVsql Sl^ x g ; 7 4 (1~P\ Io <<;L 



(3.31) 
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In conclusion, in the ERSE configuration (<fr = tt/2), the situation is in some sense the reverse of the ESR 
scheme (<p = 0). In the former the bandwidths are much larger than in either the ESR of the conventional 
interferometer. However, the more broadband curves are obtained at the cost of losing sensitivity in the 
frequency range ~ 70 — 250 Hz and this explains why the maximum signal-to-noise ratio for inspiraling 
binaries, 

S BSE ^ 1-096 for p-0.48 and I = J SQL . (3.32) 

{Of iV Jconv 

is not very different from that of a conventional interferometer. Finally let us observe that our two extreme 
cases are linked mathematically by taking p — ► — p ( JC — > K) and exchanging the two quadratures. For 
much further analysis and detail of the ERSE and ESR configurations, see Refs. HQ. 



IV. STRUCTURE OF RESONANCES AND INSTABILITIES 



We now turn our attention from the well known extreme configurations, for which previous analysis gave 
correct predictions, to the more general case < (j) < tt/2. As Figs. |, | show, the noise curves for a SR 
interferometer with frequency independent homodyne detection generically exhibit resonant features that 
vary as I , p, 4> and £ are changed. These resonances are closely related to the optical- mechanical resonances 
of the dynamical system formed by the optical field and the mirrors. A thorough study of this system must 
investigate explicitly the motion of the mirrors, instead of including it implicitly in the formulae as we did 
in this paper. It can be most clearly worked out using the formalism of linear quantum measurements 0], 
which we have recently extended to SR interferometers [ [llj . In this section, we limit our investigation to 
the resonant structures in the amplitudes of the optical fields, and for simplicity we work in the limit of a 
totally reflecting SR mirror, i.e. p = 1. This limit provides simple analytical expressions for the resonant 
frequencies as functions of the SR detuning phase </> and the light power I Q . We shall comment on the general 
case p ^ 1, which we have tackled at length in Ref. jnj, only at the end of this section. 



A. Resonances of the closed system: p — 1 

We shall investigate the free oscillation modes of the whole interferometer when the GW signal is absent 
[h(Q) = 0] and there is no output field (p — 1), so the system is closed. We consider the regime of 
classical electrodynamics, i.e. we work with the two classical quadrature fields E\ and E%, satisfying the 
same equations of motion as the quantum- field operators c\ and ci (see Fig. [l]). We shall evaluate the 
stationary modes, notably the eigenmodes and eigenvalues of the whole optico-mechanical system made of 
the end mirrors and the signal recycled optical field. We achieve this by propagating the in-going fields E\ 
and Ei (entering the beam splitter's dark port) into the conventional interferometer, along a complete round 
trip, and then through the SR cavity back to the starting point. The round-trip propagation leads to the 
following homogeneous equation for the eigenmodes: 

«20 - S in2«A 2i/3 / 1 0\ _ 1 fE x 
sin20 cos20 / \-K, l) \ \E 2 



-0, (4.1) 
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which can be simplified into the more interesting form: 

T { o e2t{ _ a+m _ 1 )T- 1 ( E J=0, 2a = arccos [cos 20 + - sin 20 j , (4.2) 

where T is a matrix whose precise form is unimportant. Note that the definition of the function arccos 
ensures that 3?(2a) ranges from to tt. The free oscillation condition is then given by: 

cos 2/3 rcs . = cos 2a — cos 2(f> + — sin 2(f> . (4.3) 



Solving Eq. (4.3) explicitly in terms of the frequency O, we obtain the rather simple analytical equation for 
the position of the resonances: 



tan 2 (p ± a / tan 4 <j> — - — — tan i 
V -*SQL 



7 2 2 

This equation is characterized by three regimes (0 < < 7r): 

• 4> > 7r/2: one real and one imaginary resonant frequency; 

• arctan[(4/ //sQL) 1 ^ 3 ] < 4> < ft/2: two real resonant frequencies; 

• < (f) < arctan[(4J //sQL) 1 ^ 3 ]: two complex conjugate resonant frequencies. 



(4.4) 



Equation (4.4) is very similar to the resonance equation that Braginsky, Gorodetsky and Khalili have derived 
for their proposal "Optical bar gravitational wave detectors" (see Appendix D of Ref. f|). 

For very low light power, I Q <C ^sql, the second term under the square root on the RHS of Eq. (4.4) goes 
to and the four roots tend to Q, = (double root) and f2 = ±7tan^. We interpret this limit as follows 
(see Ref. [0 for further details): When the coupling between the motion of the mirror and the optical 
field is zero (/„ — + 0), the resonant frequencies of the entire system are given by the resonances of the test 
mass, i.e. the free-oscillation modes of a test mass (O = 0), plus the resonances of the optical field, i.e. the 
electromagnetic modes of the entire cavity with fixed mirrors, given by Q = ±7tan0 |L3]]. When the light 
power is increased toward /sql, the coupling between the free test mass and the optical field drives the four 
resonant frequencies away from their decoupled values. By analyzing the four coupled resonant frequencies, 
we can easily identify the ones with the — (+) sign in Eq. (4.4) as remnants of the resonant frequencies of 
the free test mass (optical field). (For a more thorough discussion of these results see Ref. |y]], where we 
explicitly examine the mirror motion.) 



Let us observe that Eq. (4.3) can also be obtained as follows. By expanding the noise spectral density 



(3.5) for t — > 0, we get: 



S h (n) (-2 cos 2/3 + 2 cos2 + /C sin2 0) 2 1 



+ O(t ). (4.5) 



/ig QL (O) 8fC [cos 2 (3 (sin 2 C - cos 2 <f>) + cos 2 cos 2 C] 

The leading term of the expansion goes to zero when 2 cos 2(f)— 2 cos 2/3 rcs . + K. sin 2^ = 0, which is exactly the 
resonant condition (4.3) for the closed system derived above. This means that for (open) SR interferometers 
with highly reflecting SR mirrors, the dips in the noise curves agree with the resonances of the closed system. 
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FIG. 5. Log-log plot of ^Shity/S^i-y) versus Q/y for I D = 7sql, p = 0.95 and ( = (i.e. the second quadrature 
&2 is measured). The detuning phase (f> takes the values (going from right to left): 7r/2 — 0.19, 7r/2 — 0.39 and 

These lines 



7r/2 — 0.59. The vertical grid lines have been drawn by using Eq. (4.4) and taking the real part of fl T . 
agree well with the positions of the resonant dips. 



In practice, the real part of the resonant frequencies (4.4) for the closed system turns out to be a good 
approximation to the positions of the valleys in the noise spectral density of an (open) SR interferometer 
with high SR-mirror reflectivity. To illustrate this fact, in Fig. [| we plot the noise curves y/S^JTj for the 
second quadrature &f with I = Isql, P = 0.95 and varying <fi. The vertical lines have been drawn by solving 



Eq. ( |4.4|) numerically for f2 and taking its real part, i.e., the real part of the resonant frequencies of the 
closed systems. There is indeed very good agreement. This suggests that the gain in sensitivity comes from 



a resonant amplification effect; sec the discussion at the end of the Sec. |IV C . 

If the imaginary part of the resonant frequency is positive (negative) then, with our convention for the 
Fourier transform, the solution is unstable (stable). The best noise sensitivity curves have detuning phase 
4> in the range arctan[(4/ //sQL) 1 ^ 3 ] ^5 4>^ tt/2, which for p = 1 correspond to two real resonant frequencies, 
and no instability. However, as soon as we allow the transmissivity of the SR mirror r to be different from 
zero (as it must be in a real interferometer), we always find that one of the two resonant frequencies has a 
positive imaginary part p"l| . A more detailed analysis of the dynamics of the system has shown that this is 
a rather weak instability which typically develops on a time scale of < O.I7 and can be cured by introducing 
an appropriate control system pT| . 



B. Semiclassical interpretation of resonances for small /C: pure optical resonances 



In this section we shall focus on the optical-field resonances and shall relate our results to previous 
semiclassical analyses of SR interferometers jl3|,|l4| . 
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The test-mass motion affects the optical fields through the term JC = 2 (I //sql) 7 4 / (S1 2 (S1 2 + 7 2 )) , where 
the factor I /Isql can be considered a measure of the strength of the coupling. The quantity JC governs both 
the resonant condition and the relative magnitude of shot noise and radiation-pressure noise. In particular, 
when JC is very small, Eq. (L3) simplifies to cos 2(f) — cos2/3 res = 0, which can be solved easily, giving: 

2(±/3 ros . + <f>) = 27rn , i.e. ft Tes , = ±7 tan (f> , (4.6) 



with n an integer. Equation d4.6| ) can be explained with a simple optics argument: The quantity ±2/3 is 
the phase gained by the upper and lower GW sidebands while in an arm cavity, while <f> is the phase gained 



when traveling one way down the SR cavity. Thus 2(±/3 + cf>) is just the round-trip phase, and Eq. (4.6) 
is the resonant condition for the entire (closed) interferometer. Hence, the presence of JC in the resonant 
condition ( |4.3| ) provides the deviation from a pure optical resonance. Moreover, JC is also an indicator of 
the different scalings of I Q and m in the final expressions for the noises, and therefore it governs the relative 
magnitude of the shot noise and radiation-pressure noise - the smaller the JC, the more important the shot 
noise compared to radiation pressure noise. When JC is small, a semi-classical argument helps to explain the 
features of our noise curves. If we are close to the resonance, then feeding back the signal at that frequency 
increases the peak sensitivity while decreasing the bandwidth. Different schemes of such narrow-banding 
have been proposed, e.g., see Drever |29|. The scheme discussed here, in which the signal at the dark port is 
fed back into the arm cavities, is called signal recycling (in the narrower sense), and was invented by Meers 
p3| . If, on the other hand, we are far enough from the resonances, sideband signals are not encouraged to go 
back into the interferometer; in particular, at |/3 a nti-ros. | — |/3 r0 s. ± 7r /2|, there is antiresonance, and the signal 
is encouraged to go out. This is what is generally called resonant sideband-extraction and was invented by 
Mizuno 0, see Sec. [II C . The range in between, f3 Ics . < < /3 a nti-ros., is called 'detuned' signal recycling 
and has recently been demonstrated experimentally on the 30 m laser interferometer at Garching, Germany 
by Freise et al. |30| and at Caltech on a table-top experiment by Mason |3l|j . 

As an example of resonance (not anti-resonance), we plot in Fig. ^ the spectral density Sh(Q) when the 
second quadrature is measured, for very low light power I a — 10~ 4 /sql arid high reflectivity p — 0.95, 
and for various values of the detuning phase <f>. The vertical grid lines in Fig. |^ are drawn according to 
Eq. (L6) and indeed, there is excellent agreement. 

It is interesting to note that although for LIGO-II I — Isql, there is still a frequency band where JC is 
relatively small. This is due to the fact that JC drops very fast as SI increases. In that frequency band the 
semiclassical formalism gives a correct result for the optical resonances pS[ ], However, since the semiclassical 
approach does not take into account the motion of the arm-cavity end mirrors, it can only describe one 
resonance (and not two) in the entire spectrum. 



C. Quantum mechanical discussion of the general case: two resonances and p 7^ 1 

The correspondence between the optical-mechanical resonances and the minima of the noise curves suggests 
that the gain in sensitivity comes from a resonant amplification of the input signal, i.e. of the gravitational 
force acting on the mirrors, as already observed for optical bar GW detectors by Braginsky's group |J. Let 
us discuss this point more deeply. 
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FIG. 6. Log-log plot of ^5h(n)/S^ L (7) versus fi/7 for £ = (i.e. is measured) and for extremely low light 
power and high reflectivity: I = 10 -4 7sql and p = 0.95. <f> takes the values (going from right to left): ty/2 — 0.19, 
7r/2 — 0.39, 7r/2 — 0.59, n/2 — 0.79 and n/2 — 0.99. A series of resonances appear whose positions agree with the 



vertical grid lines drawn according to Sl rC s./7 = | tan0| [Eq. (4.6)] 



The quantum part of the input-output relation (2.2C) (with |$| <C 1 as we have assumed throughout this 
paper) reads 

M 



br * = ?_£*L ajt i,j = 1,2. (4.7) 



We find it convenient to renormalize the quantum transfer matrix: 

so det M.ij — 1. Note that this A4ij is normalized with respect to unit quantum noise. Because the CV,- are 
real, the matrix M. depends on three real parameters and we can always decompose it into two rotations 
R{9), R(<p) and a squeeze S(r) (see for details Rcf. ||), e.g., M = R(9) R(<p) S(r) R(-<p), with 

/ cos 6 - sin 6 \ ( e r \ 

R{9) = • n *)> S{r) = n -r ) ' ( 4J ) 
V sm cos J \ e / 

where the factor e r describes the stretching (r > 0) or squeezing (r < 0) of the quantum fluctuations in the 



quadrature bi [see Eqs. (4.7), (fLq)]. Note that classical optical fields always have a zero squeeze factor. 

To express the squeeze parameter r in terms of the physical parameters describing the SR interferometer, 
we simply take the trace of the matrix M. obtaining 

e- + e- = 2 + ^. (4.10) 

Hence, in a SR interferometer the squeezing (generally called ponderomotive squeezing) is induced by the 
back-action force acting on the mirror through the effective coupling /C. In particular, for small IC, we have 
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e 2r + e~ 2r ps 2 and the squeeze factor r goes to 0, which means the output field is classical. For our discussion 
below the specific expressions of 6 and <p in terms of the physical parameters are unimportant. 

From the previous discussions and the results derived in Q we have learned that the zeros of M(Q) are 
the resonant frequencies of the optical-mechanical system and the valleys of the noise spectral densities are 
their real parts. It is straightforward to show that for f2 equal to the real part of the resonances \M\ oc r 2 . 
Hence, on resonance, for typical values of the physical quantities I , p and </>, the RHS of Eq. ( 4.10 ) goes to 
a constant when r — > 0. This means that the squeeze factor r does not grow much around the resonances. 
On the other hand, the absolute value of the output signal strength [the term involving h in Eq. ( 2.20Q ], is 
given by 



V2Xt|A| h 
h S QL \M\ 

1/r , when r 



E2. 



(4.11) 



and because on resonance 1/|M 
amplification becomes stronger and stronger as r 



• the classical signal is resonantly amplified and the 
(closed system). 

This means that, by contrast with QND techniques based on static correlations between shot noise and 
radiation-pressure noise Q , (^] , in SR interferometers the ponderomotive squeezing does not seem to be the 
major factor that enables the interferometer to beat the SQL. Indeed, whereas the amplitude of the classical 
output signal is amplified near the resonances, the nonclassical behavior of the output light is not resonantly 
amplified. Therefore, the beating of the SQL in SR interferometers comes from a resonant amplification of 
the input signal: the whole system acts like an optical spring, 13 as we have described more thoroughly in 
PI, and it was also derived for optical bar GW detectors by Braginsky's group 0. 



V. INCLUSION OF LOSSES IN SIGNAL-RECYCLING INTERFEROMETERS 

In this section we shall compute how optical losses affect the noise in a SR intcferometer using the lossy 
input -output relations for a conventional interferometer [jfj and doing a similar treatment of losses in the 
SR cavity. We shall continue to use our extension of the KLMTV's formalism as developed in Sec. [0]. In 
jnj we show that when losses are included a suitable control system can be implemented to circumvent the 
instabilities. 

KLMTV § described the noise that enters the arm cavities of a conventional interferometer at the loss 
points on the mirrors in terms of a noise operator, whose state is the vacuum, with quadratures n\ and n 2 - 
The resulting lossy input-output relations read JtJ 

di = Cl e 2l!3 (l-^r) +VS e^m, (5.1) 



2 r 



d 2 = c 2 e 21 ' 3 1 - - + V£e i0 n 2 + V2K 



2 / ~ hsQL 



e 



iff 



e 



l- i (3 + e 2 ^ 



In this sense we could refer to a signal recycled interferometer as a SPRING detector, which could also stand for 
Signal Power Recycling Interferometer Gravitational detector ! 
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(5.2) 



where e = 2C/T and £ is the loss coefficient per round trip in the arm-cavity. For L1GO-1I T and C are 
expected to be T = 0.033 and £^200 x 10~ 6 , so e^O.01. The quantity £ which appears in Eqs. (5.1) and 



(5.2) is frequency dependent and is given by 

£ = 



2c 



(5.3) 



i + Wi) 2 

In the present analysis, as in Ref. (7), we do not take into account losses coming from the beam splitter. We 
expect their effect to be small compared to the losses introduced by the SR cavity and the photodetection 
process. Fig. [7] sketches the way we have incorporated losses. We describe the loss inside the SR cavity by 
the fraction of photons lost at each bounce of the interior field off the SR mirror, Asr, and we introduce 
associated noise quantum operators pi (i = 1, 2) into the inward-propagating field operator at the SR mirror 
(see left panel of Fig. 



Equations (2.18) then become 



pi 



\J\ - A S r (toi + p fx) + \J A S r Pi , e 2 = \/l - A S r (t a 2 + p f 2 ) + \/^sb,P2 • 



(5.4) 



and the noise operators pi satisfy the commutation relations ( |2.8| ) . We also assume that the state of pi is the 
vacuum. We include the losses of the photodetection process in an effective way, by modifying the output 
field operators and introducing another noise field qi with i = 1,2 (see right panel of Fig. [?]): 



b\ = y/l- A PD (r /i - p oi) + VA 



pd qi , 



A PD (r f 2 



pa 2 ) 



pd qi 



(5.5) 



Here, Apo is the photodetector loss. The noise quadrature fields qi describe additional shot noise due to 
photodetection and are assumed to satisfy Eq. fl2.8| ) and to be in the vacuum state. Following the procedure 
described in Sec. ||, we derive from Eqs. (5.1), ( |5.2| ), (5.4) and (5.5) the following input-output relations for 
the lossy SR interferometer (for simplicity we set $ = 0): 
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(5.6) 



where, to ease the notation, we have defined 



M L = 1 + p 2 e 4 ^ 



2/> j cos 20+ — sin 20 



+ A S rP -pe 



,2iP 



cos 20 



— sm20 



.-,2i£* 



-ep 



2 cos z (3(-p e 210 + cos 20) + ^ (3 + e Mt> ) sin 20 



K 



2Z0X 



(5.7) 



Note that M L , like M in Eq. ( 2.21 ), has zeros in the lower- and upper- half complex £1 plane. Hence, the 
lossy SR interferometer, like the lossless one, also suffers from instabilities. Nevertheless, we have shown in 
jllj that an appropriate control system can cure them, as in the lossless case. In the following equations we 
give the various quantities which appear in Eq. ( |5.6| ) accurate to linear order in e and Asr but to all orders 
in Apd. (We expect Asr~0.02 and Apd^O.1 |2S|| .) The various quantities read 
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FIG. 7. Sketchy view of the lossy signal-recycling interferometer. Optical losses in the signal-recycling cavity (on 
the left) are described by the noise quadratures pi, while losses due to the photodetection process (on the right) are 
included through the noise quadratures qt . 
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D\ = y/l — Apd j-(l + pe 2 ^) sin0+ie [3 + p + 2pe 4 ^ 
+e 2l/3 (1 + 5p)] sin0+iA SR e 2l/3 psin0| , 

£>2 = \/l - a pd |-(-l + pe 2 ^) cos0+ [-3 + p + 2pe 
+e 2 ^(-l + 5p)] cos0+iAs R e 2l %cos0j ; 



(5.9) 



Pn = P 2 2 = \ y/l - Apd V^sr t (-2p e 2l/3 + 2 cos 20 + K sin 20) , 
P12 = — \/l — Apd \J Asr r sin0 (2 cos + /C sin 0) , 



P<2,\ = \A — Apd a/ Asr t cos (2 sin — /C cos < 
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(5.12) 



Similarly to Sec. Ill A, we follow KLMTV's method to derive the noise spectral density of a lossy SR 



interferometer [see Eq. (3.5)] 
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Exploring numerically this equation, we find that for the loss levels expected in LIGO-II (e~0.01, Apd~ 
0.1, Asr ~ 0.02 |2l[]), the optical losses have only a modest influence on the noise curves of a lossless SR 
interferometer. For example, in Fig. || we compare the lossless noise spectral densities with the lossy ones 
for the two quadratures b\ and 6 2 . The main effect of the loss is to smooth out the deep resonant valleys. 
More specifically, for (i) the physical parameters used in Fig. [|, (ii) a net fractional photon loss of 1% in the 
arm cavities (e = 0.01) and 2% in each round trip in the SR cavity (Asr = 0.02), and (iii) a photodetector 
efficiency of 90% (Apd = 0.1), we find that the losses produce a fractional loss in signal-to-noise ratio for 
inspiraling binaries [see Eqs. (3.8), (|3.9|)] of 8% and 21%, for the first and second quadratures, respectively. 



The reason why we get a modest effect from optical losses as compared to schemes using squeezing or FD 
homodyne detections^ rests on the fact that our gain in sensitivity mostly comes from resonant amplification, 
which is much less susceptible to losses than quantum correlations. This general consideration has long been 
understood by Braginsky, Khalili and colleagues and underlies their motivation for the GW "optical bar" 
detectors |)). 

VI. CONCLUSIONS 

In this paper we have extended the quantum formalism recently developed Q| for conventional interfer- 
ometers (LIGO-I/TAMA/Virgo), to SR interferometer such as LIGO-II. The introduction of the SR cavity 



14 Note that in KLMTV |?J they assumed a loss factor for end-mirrors which is 10% of our value, and they also did 
not take into account losses coming from the photodetection. 
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has been planned as an important tool to reshape the noise curves, making the interferometer work either 
in broadband or in narrowband configurations. This flexibility is expected to improve the observation of 
specific GW sources |L7) . Quite remarkably, our quantum mechanical analysis has revealed other significant 
features of the SR cavity. 

First, the SR mirror produces dynamical correlations between quantum shot noise and radiation-pressure- 
fluctuation noise which break the light's ability to enforce the SQL of a free mass, allowing the noise curves 
to go below the SQL by modest amounts: roughly a factor two over a bandwidth A/~/. Before our work, 
researchers were unaware of the shot-noise / radiation-pressure correlations and thus omitted them in their 
semiclassical analysis of the straw-man design of LIGO-II . The goal of beating the SQL in LIGO-II can 
be achieved only if all sources of thermal noise can also be pushed below the SQL and indeed much R&D will 
go into trying to push them downward. It turns out that even with current estimates of the LIGO-II thermal 
noise pql , which are a little above the SQL, the net noise (thermal plus optical) is significantly affected by 
the shot-noise/ radiation- pressure correlations. Indeed, the correlations lift the noise at low frequencies 
10 Hz < fi/27r < 50 Hz, as compared to the semiclassical estimations, even though in this frequency range the 
optical noise may already be very much larger than the SQL. This is due to the inaccuracy of the semiclassical 
method in estimating the effect of the radiation-pressure force, which is important in this region. In the 
middle frequency range, i.e. near 100 Hz, the SQL-beating effect cannot lower the total noise much because 
of the thermal contribution. The effect of the correlations in the implementation of LIGO-II will be clarified 
and sharpened once the readout scheme has be specified [ pi[ . 

Second, we have learned that the dynamical correlations arise naturally from the nontrivial coupling 
between the antisymmetric mode of motion of the four arm-cavity mirrors and the signal recycled optical 
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fields. This dynamical coupling invalidates the naive picture, according to which the arm cavity mirrors are 
subject only to random quantum-vacuum fluctuations. The SR interferometer responds to a GW signal as 
an optical spring jll) , and this oscillatory response gives the possibility for resonant amplification of the GW 
signal. The optical-mechanical system is characterized by two resonances and one of them is always unstable, 
so a control system must be introduced to stabilize it [fill . In the limit of a highly reflecting SR mirror we 
have worked out analytically a very simple equation which locates the positions of the resonant frequencies. 
Whereas the amplitude of the classical output signal is amplified near the resonances, the quantum noise 
is not particularly affected by them. All this suggests that the beating of the SQL in SR interferometers 
comes primarily from the resonant amplification of the input GW signal, as also occurs in "optical bar" GW 
detectors (9). 

The inclusion of losses does not greatly affect the SR interferometer. This is due to the fact that the 
improvement in the noise curves rests primarily on a resonant amplification and only modestly on pon- 
deromotive squeezing. It is worthwhile to point out that the SR interferometers bears strong similarity to 
the "optical bar" detectors proposed by Braginsky, Khalili and colleagues Both of them can be viewed 
as oscillators with two different eigenfrequencies. However, because in SR interferometers the light plays 
the double role of providing the restoring force and being a probe to monitor the mirror displacements, 
we are forced to introduce in SR interferometers much higher laser power, to circulate in the arm cavi- 
ties (~ lMWatt), than in the "optical bar" scheme. Nevertheless, like the "optical bar" scheme, the SR 
interferometer is still less susceptible to optical losses than many other schemes designed to beat the SQL. 

It is now important to identify the best SR configuration, i.e. the choice of the physical parameters (light 
power I Q , SR detuning <f>, reflectivity of SR mirror p, quadrature phase £, and the read-out scheme: homodyne 
or modulation/demodulation) that optimizes the signal-to-noise ratio for inspiraling binaries, for low-mass 
X-ray binaries, and for other astrophysical GW sources. We shall discuss this issue in a forthcoming paper 
II- 

Finally, our analysis has shown that dynamical correlations, i.e. correlations that are intrinsic to the 
dynamics of the test mass-optical field system (i.e. they are not due to specific read-out schemes, as in the 
case of homodyne detection on a conventional interferometer), are present when the carrier frequency uj is 
detuned from resonance (<f> ^ 0) or antiresonance (0 ^ tt/2) in the SR cavity. This suggests a speculation 
that it could be worthwhile to investigate a LIGO-II configuration (see Table |) without a signal recycling 
mirror, in which the correlations are produced by detuning the arm cavities. However, this case will require 
a very careful analysis of the radiation-pressure force acting on the arm-cavity mirrors [ p0[ , JlS| ] . 
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APPENDIX A: REMARK ON COMMUTATION RELATIONS AMONG QUADRATURE FIELDS 
IN CAVES-SCHUMAKER TWO-PHOTON FORMALISM 



As originally pointed out by Braginsky's group Q and discussed by BGKMTV |16|, the output variables 
of the GW interferometer should commute with themselves at different times, to guarantee that no other 
quantum noise is necessarily introduced into the measurement result once further manipulations are per- 
formed on the output. Indicating generically by 0(t) the output quantity, the following conditions should 
be satisfied, 



p(t),o(t')] = o vt,i', [0(Q),o Jf (n')] = vn,n'. 



(Al) 



If we assume that the system's output is one quadrature of the quantized electromagnetic field (EM) [see 
Eq. ( 2.1C| )1, with the GW signal encoded at side-band frequency Jl around the carrier frequency ui , then 
the presence of terms proportional to £1/lu in Eq. (2.5) prevents the output quadratures from commuting 
with themselves at different times. However, Braginsky et al. [[u| anticipated that, in the case of LIGO- 
I/TAMA/Virgo, the quadrature fields at the dark port should anyway satisfy very accurately the Fourier- 
domain condition given by Eq. (Al), because the side-band frequency f2 (1Hz < £l/2ir < 10 3 Hz) is much 
smaller than the carrier frequency oj {uj ~ 10 s _1 ). In this Appendix we investigate this approximation 
in much more detail, estimating the amount of extra noise which will be present in the final noise spectral 



density as a result of condition (Al) being violated. Henceforth, for simplicity we restrict our analysis to 
conventional interferometers. 

If the readout scheme is implemented by photodetection, then only a small frequency band around lu con- 
tains the final output signal. Hence, it is physically justified to introduce a cutoff A in the frequency domain 
which automatically discards all the Fourier components of the EM field outside the range [uj — A, oj + A] 
with < A < u . As a consequence, Eq. ( f2~6| ) for the EM field can be rewritten as [see also Eqs. (4.22) of 
Ref. ©] 
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where a+(fi) and a_(Sl), with SI < A, are defined by Eq. ( |2.3| ) and the rescaled quadrature fields Of{t) are 
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with the quadrature operators given by Eq. ( |2.7D . Evaluating the commutation relations among the quadra- 
ture operators we find [see also Eqs. (4.31) of Ref. [p2[]: 
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Note that Eq. (A5) differs from the one appearing in Eq. ( |2.S| ), where we approximated dj and a\, as 
commuting. The non-vanishing commutation reiations in Eq. ( |A5|) explicitly yield a non-vanishing two-time 
commutator for Of". In particular, a straightforward calculation gives (i = 1,2): 



cW(t,f)= [oftf),0ftO] =i— 

1 1 w 



At cos(At) — sin(Ar) 
7T (At) 2 



T = t-t' . 



(A7) 



Therefore Of(t) cannot be the final output and there must be some unavoidable additional quantum noise due 
to the fact that Of(t) has a non-vanishing two-time commutator. In LIGO-I/TAMA/ Virgo this additional 
noise is introduced in the output during the final process of photodetection. A more detailed study would 
involve a very technical analysis of the photodetection's dynamics, but fortunately, as we shall see in the 
following, a simple estimation of the order of magnitude of this additional quantum noise suggests that it is 
very small and we can realistically neglect it. 

We find it convenient to estimate the additional quantum noise by calculating the noise induced by the 
photodetector approximated as a linear measurement device coupled to the quadrature fields. Having 
fixed the cutoff frequency A and working in the Fourier domain, we can write the final output as: 



orm = om + z^m + R ^m prw 
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The last two terms in Eq. ( A8) are the shot noise and the back-action noise of the photodetector (PD) and 
describe the efficiency and the strength of perturbation of the PD on the quadrature field, respectively. Let 
us assume that there is no correlation between Zf B and Ff u . Hence, Zf D and Ff T> satisfy the uncorrelated 
version (3.16) of the uncertainty relation, that is 



S gPIi gPD Sp PDjfPD > H . 



(A10) 



We are interested in evaluating the overall quantum noise. We first write the output in the form Signal 
Noise as 



(All) 



where Vih is the part of Of~(fl) that contains the signal, while Qf (£1) contains the quantum fluctuations. 
Using Eq. (All), the overall noise spectral density is (i = 1,2): 
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Here we are assuming that as a consequence of the homodyne detection, the EM field impinging on the pho- 
todetector is composed of carrier light plus quantum fluctuations, and thus the light intensity measured by the 
photodetector is linear in the annihilation and creation operators. 
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The first term in Eq. ( |A12| ) describes the quantum noise of an interferometer when the non-vanishing 
commutators of the quadrature fields have been ignored and ideal photodetection is applied. The second 
term in Eq. ( A12 ) describes the additional shot noise introduced by the photodetection process. Finally, the 
third term comes from the back-action force acting on the measured quadrature (i = 1 or 2) because it does 
not commute with itself at different times. Let us notice that, given Eq. ( A10 ), the second and third noise 
contributions appearing on the RHS of Eq. ( A12 ) are complementary. Indeed, the larger the shot noise, 
the weaker the minimum force the photodetector must apply to the quadrature fields and the smaller the 



back-action noise. More specifically, there is a lowest achievable value for the PD part in Eq. (A12) given 
by: 
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Using Eq. (gj) we derive 1/|S;| 2 = h 2 sc(L /2K and S q a q a = [K 2 + 1) > 1. Recalling that 10Hz < SI/2tt < 
10 3 Hz and ui ~ 10 15 sec -1 , fixing A to a value larger than the typical St, e.g., A/2ir^ 10 MHz, and adjusting 
the PD such that S z pt> z pt> and S f p-d f p-d satisfy the minimal uncertainty relation [the equality sign in 
Eq. ( A10 )1, we find that the minimal achievable PD noise is ~ 1CP 7 times the conventional shot noise. 
Therefore, we can totally ignore the quantum noise introduced by the fact that the quadrature fields do not 
commute with themselves at different times in Eq. (|A5|). Note the importance of the cutoff A. If we had 
taken A~ oj , the limit on the PD noise would have been of the same order of magnitude as the shot noise for 
a conventional interferometer and it would not have been realistic to neglect the quantum noise introduced 
by the non- vanishing commutation relations of the quadrature fields. 

So far we evaluated the minimum quantum noise that the photodetector, coupled linearly to the quadrature 
field, can introduce. Let us now try to give a realistic value of it. To estimate S z po z pd, we can just recall 
that in the case of a lossy photodetector we have (see the discussion of lossy interferometers in Sec. ^) 



PD Pi 



(A14) 



where pi with i = 1,2 are quadrature operators in the vacuum state. We expect Apd ~ 0.1 — 0.2, hence 



?d > 10 -2 x S. 



shot noise 



, which is five orders of magnitude larger than the lowest achievable limit 



discussed above with A/27r~10MHz. Therefore, if one can justify fixing the cutoff A/27T at 10 MHz, and if 
the uncertainty relation ( |A10| ) is satisfied with the equality sign, then one can conclude that the inefficiency 
will dominate over the minimum possible back-action noise by five orders of magnitude. Hence, we are 
justified in disregarding the non- vanishing two-time commutators of the quadrature fields in Eq. (A5). 
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